Molecular scale simulations on thermoset polymers: A review

This article reviews the field of molecular simulations of thermoset polymers. This class of polymers is of interest in applications ranging from structural components for aerospace to electronics packaging and predictive simulations of their response is playing an increasing role in understanding the molecular origin of their properties and complementing experiments in the search for tailored materials for specific applications. It focuses on modeling and simulation of the process of curing to predict the molecular structure of these polymers and their thermomechanical response by all-atom molecular dynamics simulations. Results from Monte Carlo and coarse-grained simulations are briefly summarized. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2015 , 53, 103–122     

Free volume characteristics of 2,2-bistrifluoromethyl-4,5-difluoro-1,3-dioxole-co-tetrafluoroethylene copolymers: Effect of composition and molecular weight

2,2-bistrifluoromethyl-4,5-difluoro-1,3-dioxole-co-tetrafluoroethylene (PDD-TFE) copolymer is a good candidate to prepare gas separation membranes with excellent permeability due to its free volume characteristics. However, the influence of PDD-TFE copolymer structure on its free volume characteristics is less studied. In this paper, PDD-TFE copolymers with different compositions and molecular weights were synthesized, and their free volume characteristics were analyzed by positron annihilation lifetime spectroscopy and a molecular dynamics simulation. It indicated that the molar fraction of PDD in copolymers had a significant effect on free volume characteristics, while the molecular weight of copolymers exerted a slight influence on free volume when the molecular weight exceeded a critical region (intrinsic viscosity [η] > 68 ml g⁻¹). PDD-TFE copolymers with greater PDD molar fractions (i.e., 72% and 84%) showed bimodal distributions in positron lifetime and free volume size distributions, while PDD-TFE copolymers with lower PDD molar fractions (i.e., 27% and 35%) exhibited a single peak. The long-lifetime parameter τ₃ was assigned to micro-cavities formed by [-(TFE)_y-PDD-] segments and τ₄ was attributed to micro-cavities formed by [-(PDD)_x-TFE-] segments. The cis and trans transitions of PDD led to a local multilayer spiral structure with a 2.6–4.3 Å layer spacing, which would also increase the free volume of copolymers.     

Modification of an implicit approach based on nonstandard rules for structural dynamics analysis

The present study aims to modify a recently suggested implicit approach consisted of the approximate Euler method and closed-form exponential mapping (herein referred to as the Liu scheme) for the dynamic analysis of structures. Such modification has been developed based upon nonstandard rules. The equation of motion is formulated in the augmented dynamic space to apply the exponential mapping as a group preserving scheme. The formulation of the proposed method involves the hyperbolic sine and cosine functions. The method is therefore prone to divergence due to the behavior of the hyperbolic functions in structures with a high ratio of stiffness to mass. In the present study, to consider the properties of the structural equation into the formulation of the time step size and thereby avoid the divergence, a parameter, known as stability parameter, is thus derived from the exact solution of the equation of motion based on nonstandard rules. Embedding this parameter into the proposed method improves its stability. Afterward, for evaluating the performance of the proposed method, it is applied to several structures with different loading patterns while implemented in programing environment of the Matlab software. The results are compared to those of several commonly used numerical methods in structural applications. It is found that the proposed method has acceptable convergence and accuracy, and low time consumption compared to several commonly used methods. Furthermore, its stability is guaranteed by embedding the stability parameter into the proposed method.     

Self-similar asymptotic behavior for the solutions of a linear coagulation equation

In this paper we consider the long-time asymptotics of a linear version of the Smoluchowski equation which describes the evolution of a tagged particle moving in a random distribution of fixed particles. The volumes v of these particles are independently distributed according to a probability distribution which decays asymptotically as a power law $v^{?\sigma }$. The validity of the equation has been rigorously proved in [22] taking as a starting point a particle model and for values of the exponent $\sigma >3$, but the model can be expected to be valid, on heuristic grounds, for $\sigma >\frac{5}{3}$. The resulting equation is a non-local linear degenerate parabolic equation. The solutions of this equation display a rich structure of different asymptotic behaviors according to the different values of the exponent σ. Here we show that for $\frac{5}{3}<\sigma <2$ the linear Smoluchowski equation is well-posed and that there exists a unique self-similar profile which is asymptotically stable.     

Molecular dynamics of the intramolecular 1, 3-dipolar ene reaction of a nitrile oxide and an alkene: non-statistical behavior of a reaction involving a diradical intermediate*

ABSTRACT

Potential energy surfaces and molecular dynamics of the intramolecular 1, 3-dipolar cycloaddition and ene reaction of a nitrile oxide with an alkene were performed in the gas phase and in dichloromethane with density functional theory. One hundred trajectories were propagated in the gas phase and in dichloromethane, respectively. Twenty percent of the trajectories in the gas phase involve bicyclic intermediate and the mean time gap is 472fs. A dynamically stepwise reaction is observed. In dichloromethane, more reactive trajectories were obtained and the time gap is larger than that in the gas phase.     

Conformational Dynamics Underlies Different Functions of Human KDM7 Histone Demethylases

The human KDM7 subfamily histone H3 Nϵ-methyl lysine demethylases PHF8 (KDM7B) and KIAA1718 (KDM7A) have different substrate selectivities and are linked to genetic diseases and cancer. We describe experimentally based computational studies revealing that flexibility of the region linking the PHD finger and JmjC domains in PHF8 and KIAA1718 regulates interdomain interactions, the nature of correlated motions, and ultimately H3 binding and demethylation site selectivity. F279S an X-linked mental retardation mutation in PHF8 is involved in correlated motions with the iron ligands and second sphere residues. The calculations reveal key roles of a flexible protein environment in productive formation of enzyme-substrate complexes and suggest targeting the flexible KDM7 linker region is of interest from a medicinal chemistry perspective.     

Correlating Charge-Carrier Dynamics with Efficiency in Quantum-Dot Solar Cells: Can Excitonics Lead to Highly Efficient Devices?

The photovoltaic performance of quantum-dot solar cells strongly depends on the charge-carrier relaxation and recombination processes, which need to be modulated in a favorable way to obtain maximum efficiency. Recently, significant efforts have been devoted to investigate the carrier dynamics of nanocrystal sensitizers, both in solution and deposited on TiO₂ photoanodes, with the aim to correlate the excitonics with solar-energy conversion efficiency. This Minireview summarizes some proof of the concepts that efficiency can be directly correlated to the exciton dynamics of quantum-dot solar cells. The presented findings are based on CdSeS alloy, CdSe/CdS core/shell, Au/CdSe nanohybrids, and Mn-doped CdZnSSe nanocrystals, where the favourable excitonic processes are optimized to enhance the efficiency. Future prospects and limitations are addressed as well.     

Mass distribution and skewness for passive scalar transport in pipes with polygonal and smooth cross sections

We extend our previous results characterizing the loading properties of a diffusing passive scalar advected by a laminar shear flow in ducts and channels to more general cross‐sectional shapes, including regular polygons and smoothed corner ducts originating from deformations of ellipses. For the case of the triangle and localized, cross‐wise uniform initial distributions, short‐time skewness is calculated exactly to be positive, while long‐time asymptotics shows it to be negative. Monte Carlo simulations confirm these predictions, and document the timescale for sign change. The equilateral triangle appears to be the only regular polygon with this property—all others possess positive skewness at all times. Alternatively, closed‐form flow solutions can be constructed for smooth deformations of ellipses, and illustrate how both nonzero short‐time skewness and the possibility of multiple sign switching in time is unrelated to domain corners. Exact conditions relating the median and the skewness to the mean are developed which guarantee when the sign for the skewness implies front (more mass to the right of the mean) or back (more mass to the left of the mean) “loading” properties of the evolving tracer distribution along the pipe. Short‐ and long‐time asymptotics confirm this condition, and Monte Carlo simulations verify this at all times. The simulations are also used to examine the role of corners and boundaries on the distribution for short‐time evolution of point source , as opposed to cross‐wise uniform, initial data.     

Studying the Electrically Driven Switching of the Planar Light Guide

Liquid crystalline materials can be used as an active media for the new generation of planar light guides. The main characteristics which governs light-guiding and switching abilities of such devices are the spatial distributions of the refraction indices (defined via the distribution of nematic director) for the liquid crystal confined within a light-guiding pore. We aim to obtain these distributions from the molecular dynamics simulation of the liquid crystalline cell with the homeotropic boundary conditions being applied. We discuss the reorientation kinetics of the homeotropic-planar transition and obtain the equilibrium director profile upon application of the planar reorienting field.